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We provide a thorough exposition, including technical and numerical details, of previously 
published results on the quantum stabilization of cosmic strings. Stabilization occurs through 
the coupling to a heavy fermion doublet in a reduced version of the standard model. We 
combine the vacuum polarization energy of fermion zero-point fluctuations and the binding 
energy of occupied energy levels, which are of the same order in a semi-classical expansion. 



' Populating these bound states assigns a charge to the string. We show that strings carrying 



fermion charge become stable if the electro-weak bosons are coupled to a fermion that is less 
^ ' than twice as heavy as the top quark. The vacuum remains stable in our model, because 

neutral strings are not energetically favored. These findings suggests that extraordinarily 
large fermion masses or unrealistic couplings are not required to bind a cosmic string in the 
standard model. 



I. INTRODUCTION 

O ■ It is well-known that the electroweak standard model and many of its extensions have the 

potential to support string-like configurations that are the particle physics analogs of vortices or 
magnetic flux tubes in condensed matter physics. Such objects are usually called cosmic strings 
to distinguish them from the fundamental variables in string theory, and also to indicate that they 
' typically stretch over cosmic length scales. 

The topology of string-like configurations is described by the first homotopy group 
where ^ is the manifold of vacuum field configurations far away from the string. In typical 
electroweak-like models, a Higgs condensate breaks an initial gauge group G down to some sub- 
. group H, so that ^ ^ G/H. Topologically stable strings are therefore ruled out in the electroweak 

standard model SU{2) x C/(l) — U{1) because G/H is simply connected. Nevertheless, one could 
envision a GUT and/or supersymmetric extension in which a simply connected group G breaks 
down to the electroweak SU{2) x U{1) at a much higher scale, so that Ili{G/{SU{2) x U{1))) 
^ I is nontrivial. Since such GUT strings would have enormous energy densities, they could be seen 

by direct observation using gravitational lensing [l|, 0] or by signatures in the cosmic microwave 
background [1]. Moreover, a network of such strings is a candidate for the dark energy required to 
explain the recently observed cosmic acceleration [3, . 

The absence of topological stability does not imply that electroweak strings (or Z-strings 0-0] ) 
are unstable or irrelevant for particle physics. While their direct gravitational effects are negligible, 
Z-strings can still be relevant for cosmology at a sub-dominant level [^,E3]- Their most interesting 
consequences originate, however, from their coupling to the standard model fields. Z-strings 
provide a source for primordial magnetic fields [8| and they also offer a scenario for baryogenesis 



with a second order phase transition Uj, In contrast, a strong first order transition as required 



by the usual bubble nucleation scenario is unlikely in the electroweak standard model [ij] , without 



non-standard additions such as supersymmetry or higher-dimensional operators . Because the 
core of the Z-string is characterized by a suppressed Higgs condensate, it allows for both the 
copious baryon number violation and the out-of equilibrium regions required by the Sakharov 
conditions, without relying on a first order phase transition. 

However, these interesting effects are only viable if Z-strings are energetically stabilized by their 
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coupling to the remaining quantum fields. The most important contributions are expected to come 
from (heavy) fermions, since their quantum energy dominates in the limit Nq — ?• oo, where Nc is 
the number of QCD colors or other internal degrees of freedom. The Dirac spectrum in typical 
string backgrounds is deformed to contain either an exact or near zero mode, so that fermions 
can substantially lower their energy by binding to the string. This binding effect can overcome 
the classical energy required to form the string background. However, the remaining spectrum of 
modes is also deformed and for consistency its contribution (the vacuum polarization energy) must 
be taken into account as well. Heavier fermions are expected to provide more binding since the 
energy gain per fermion charge is higher; a similar conclusion can also be obtained from decoupling 
arguments 15|, ll6|]. Dynamical stability of Z-strings in the full standard model also would suggest 
that they are presently observable. 

A number of previous studies have investigated quantum properties of string configurations. 
Naculich 17[ has shown that in the limit of weak coupling, fermion fluctuations destabilize the 
string. The quantum properties of Z-strings have also been connected to non-perturbative anoma- 
lies The emergence or absence of exact neutrino zero modes in a Z-string background and 
the possible consequences for the string topology were investigated in [191]. A first attempt at a 
full calculation of the fermionic quantum corrections to the Z-string energy was carried out in 
ref. [i^]. Those authors were only able to compare the energies of two string configurations, rather 
than comparing a single string configuration to the vacuum because of limitations arising from the 
non-trivial behavior at spatial infinity which we discuss below. The fermionic vacuum polarization 
energy of the Abelian Nielsen-Olesen vortex 21 1 has been estimated in ref. [i^] with regularization 
limited to the subtraction of the divergences in the heat-kernel expansion. Quantum energies of 
bosonic fluctuations in string backgrounds were calculated in ref. 23j] . Finally, the dynamical fields 
coupled to the string can also result in (Abelian or non- Abelian) currents running along the string's 
core. The time evolution of such structured strings was studied in ref. 2J], where the current was 
induced by the coupling to an extra scalar field. 

We have previously pursued the idea of stabilizing cosmic strings by populating fermionic bound 
states in a 2 + 1 dimensional mode l 12511 . Many such bound states emerge including, in some 
configurations, an exact zero-mode [13]. Nonetheless, stable configurations were only obtained 
for extreme values of the model parameters. In 3 + 1 dimensions, stability is potentially easier to 
achieve because quantization of the momentum parallel to the symmetry axis supplies an additional 
multiplicity of bound states. 

In this paper, we will employ the phase-shift approach, or spectral method, to compute the 
complete 0{h) fermionic contribution to the string energy from first principles. This is not a simple 
task, since the string has a vortex structure that introduces non-trivial field winding at spatial 
infinity. The standard spectral methods are thus not directly applicable since scattering theory off 
the string is ill-defined. More precisely, the Born expansion to the vacuum polarization energy, 
which in the phase shift approach is identified with the Feynman series, does not exist for the string 
background in its standard formulation. Recently we have shown how to overcome these problems 
by choosing a particular set of gauges [26], |27|. Numerical results of the full calculation of the 
string's quantum energy were first reported in ref. [2^. Here we will present the technical details 
of our calculation along with improved numerical data and a discussion of possible consequences 
of our finding. 

This paper is organized as follows: In the next section we describe our model and the string 
configuration. We then discuss the fermion Hamiltonian of our model and, in particular, how a 
local gauge transformation can be used to solve the technical problem of long-ranged gauge po- 
tentials in the string background. We also present the grand spin decomposition of the scattering 
problem. Section IIVI gives a detailed account of our method for computing the fermion quantum 



energy, which is based on the spectral approach [29l] and the interface formalism [30|]. The indi- 
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vidual contributions to the string's quantum energy are described in separate subsections, while 
some (lengthy) numerical details are deferred to appendices El [Bl and O The omission of boson 
fluctuations causes the model not to be asymptotically free which then introduces an unphysical 
Landau pole. In appendix D we verify that our results for the vacuum polarization energy are 
not affected by this artifact. In section |Vl we explain our variational search for a stable string 
configuration. To occupy fermion levels in the string background, we introduce a quantity similar 
to the chemical potential in statistical mechanics, which allows us to compute the total binding 
energy of the string as a function of the prescribed fermion charge of the string. 

Our numerical results are presented in detail within section IVII We show the parameter de- 
pendence of the individual contributions to the string's quantum energy. The stable configuration 
is discussed in more detail and it is shown that for the most stable configuration the gauge field 
contribution is negligible compared to the deformation of the Higgs field. Stabilization occurs for 
otherwise realistic parameters if the Yukawa coupling is increased by about 70% from the value 
for the top quark. Since we keep all other parameters as suggested by the standard model, this 
corresponds to a fermion mass of about 300 GeV. We close in section IVlII with a brief summary of 
our results, a discussion of its implication for the electroweak standard model and an outlook on 
possible directions for future work. 

We have published some of the results earlier [2^ and therefore focus on the technical aspects 
of the calculation here. 



II. THE MODEL 

We consider a left-handed SU (2) gauge theory in which a fermion doublet ^ = is 

coupled to a triplet gauge field Wn = 17 and a Higgs doublet (t> = { . ]■ Both 

components, and ^b, are Dirac four-spinors. This model is intended to represent the electroweak 
interactions, where we introduce some technical modifications to simplify the analysis: 

1. we set the Weinberg angle to zero so that electromagnetism decouples and the gauge bosons 
become degenerate in mass; 

2. we neglect QCD interactions, although the color degeneracy, Nq = 3, is included in the 
quantum energy arising from the fermions; 

3. we only consider a single fermion doublet and neglect inter-species (CKM) mixing and mass 
splitting within the doublet. 

With these adjustments, the bosonic part of our model is described by the Lagrangian 

£-<p,w = -^tr (G'^'^G^,) + hT{D^<P)^D,<^ - ^tr (c^tcK - , (1) 
where the Higgs doublet is written using the usual matrix representation 

$ = 



The gauge coupling constant g enters through the covariant derivative Z?^ = — i gW^ , and the 
SU{2) field strength tensor is 

Gf,, = d^W,-d,W^-ig[W^,W,]. (2) 
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We treat the bosonic fields as a classical background, ignoring the effects of bosonic fluctuations. 
This approach can be justified formally by the limit of a large number of colors Nc — oo, even 
though no QCD interactions are included: Since the quarks carry a color quantum number in the 
fundamental representation of the color group SU{Nc), their contribution to the quantum energy 
is enhanced by a factor Nc as compared to the bosonic quantum contribution. Hence we compute 
the leading quantum corrections to the classical background energy from the fermion Lagrangian 



= {PlIP +PRp)^-fi' (<PPr + <I>^Pl) ^ ■ 



(3) 



Here, Pr^l = 5 (1=1=75) are projection operators on left /right-handed components, respectively, 
and the strength of the Higgs- fermion interaction is parameterized by the Yukawa coupling /, 
which gives rise to the fermion mass, m = fv, once the Higgs acquires a vacuum expectation 



value {vev) v, where (det(<I>)) 



^ 0. All other masses in this model are also a result of the 



symmetry breaking Higgs condensate, viz. the gauge boson mass Myy = gv/^/2 and the Higgs 
mass rriH = 2v ^/X. This similarity with the standard model of particle physics suggests the model 
parameters 



9 = 0.72, 



177 GeV, mH = 140GeV, / = 0.99, 



(4) 



by taking the fermion doublet to have the mass of the top quark. Finally, the counterterm La- 
grangian necessary to renormalize the quantum energy will be listed with the computational details 
in eq. below. 

As mentioned earlier, we are particularly interested in the Z-string background configuration. 
If we consider a single straight (infinitely extended) string along the z-axis, the corresponding 
boson fields depend only on the planar polar coordinates, i.e. the distance p from the symmetry 
axis and the corresponding azimuthal angle ip. In Weyl gauge Wq = 0, we have 



■ u ^ foip) - I icos(ei)e 

W = nsm(4i) 

P V-zcos(6)e^"^ 




(5) 



/sin(6)e— ^ -icos(6)' 
^ = ^/^(^)( ..... ... 



— icos(^i) sin(^i)e*' 



(6) 



The Z~boson component = of this configuration has the familiar shape of an Abelian 
Nielsen-Olesen string of winding number n, although the entire non-Abelian configuration is 
smoothly deformable into the vacuum and thus not stable for any topological reason. We have left 
the analog of winding number n general, although we will only consider n = 1 in our numerical 
treatment below. The additional variational parameter ^1 G [0, '^/2\ was introduced to include a 
non-trivial gauge field in the string background; the same parameter also determines the orienta- 
tion of the Higgs field on the chiral circle. Then the classical energy per unit length of the string 
is a functional of the profile functions /g(p) and fnip), 



Ed 



2tt j pdp In sin ^1 



9^ 
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+ 



f2 
pp2 



(1 - fcY 



I £ff I Ph 
P 4/2 



(7) 



where the radial integration variable is related to the physical radius by Pphys = p/m and pu ^ 
mH/m. The radial functions /g(p) and fnip) in the string configuration, eqs. ([5]) and ([UD approach 
unity at large distances and vanish at the string core p = 0. Typically, they will have similar 
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FIG. 1: (Color online) Comparison between the Nielsen-Olesen profiles (full lines) and a fit using the 
variational ansatz of eq. (|S]) (dashed lines). The variational parameters are wr — 1-64/ and wg — 3.85/. 



shapes to the familiar Nielsen-Olesen string, with both W and <I> going as ^(/j) at p 
ambiguities from an undefined azimuthal angle tp. We choose a convenient form, 



fnip) = 1 - exp 



P 

WH 



and /g(p) = 1 - exp 



\wgJ 



to avoid 



(8) 



with two width parameters, wh and wg, which we also measure in inverse multiples of the fermion 
mass m. Together with the angle describing the gauge field admixture in the string, we thus have 
three variational ansatz parameters, {wh,wg,S,i), in addition to the model parameters v (which 
sets the overall scale) and the three couplings f,g and A, that are discussed above. 

To assess the quality of the variational ansatz, eq. dH) we see how well it is capable of fitting 
the Nielsen-Olesen profiles which minimize the classical energy, eq. ^ for = '^/2. As seen from 
figure [1] there is a minor discrepancy at large distances for the gauge field profile /g due to the 
Gaufiian decaying faster than any exponential function. This discrepancy affects the result for the 
classical energy in a negligible manner. For fixed = the true minimum is at 7.56f^ while the 
variational profiles yield 7.72v'^. 



III. DIRAC HAMILTONIAN 

The fermionic quantum corrections to the string background are computed in several steps. 
First, we extract the Dirac Hamiltonian associated with the Lagrangian eq. ^ and observe that 
the ansdtze, eqs. ([5]) and ([6]), do not depend on the z-coordinate (along the string symmetry axis). 
Hence this coordinate does not appear explicitly in the Hamiltonian and the z-dependence of the 
corresponding wave functions is simply e~^^''^. To compute the vacuum energy with such a trivial 



coordinate, we use the interface formalism [301 1 . which gives the quantum energy per unit length in 
terms of the two-dimensional spectrum in the plane perpendicular to the string. This formulation 
accounts for the integration over the longitudinal momentum pz using sum rules for the scattering 



data |3ll . l32l | to cope with the associated ultra-violet divergences. It then remains to solve the 
scattering problem for the Hamiltonian in the plane perpendicular to the string. 

Although we are thus left with a seemingly well-defined two-dimensional Dirac problem, the 
spectral method cannot be readily applied to compute the vacuum energy, because the long range 
of the string gauge field prevents us from setting up a well-defined scattering problem. There are 
two ways to circumvent this problem: As motivated by the study of quantum effects for QED flux 



tubes [33(], a return string was introduced in ref. [27| to unwind the gauge field at a large distance 
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from the string core. The assumption was that the energy of the return string is small when the 
unwinding is done smoothly enough and, in particular, that the associated energy density can be 
well separated from the proper string core contribution. Although these assumptions could be 
verified, the necessity to repeat the (expensive) calculation of the vacuum energy with varying 
return string positions to identify the core contribution made the return string method inefficient 
for actual calculations. 

An easier method was devised in ref. [2^ . It is based on the simple observation that the Dirac 
spectrum is gauge invariant, i.e. a local isospin rotation can be used to unwind the string gauge 
field at spatial infinity at the price of strong singularities in the origin {singular gauge), or to make 
the gauge field regular at the origin at the price of long-ranged fields at spatial infinity {regular 
gauge). The solution is to combine the good features of both gauges by means of a local gauge 
rotation that looks singular for large distances and regular for small distances. Thus, we make a 
local gauge rotation on the Dirac Hamiltonian H ^ H U with 



U = -Plcxp (i(,{p) n ■ + Pr with 



n 



cos{n(p) 
— sin(n(/7) 




(9) 



Here ^(p) is an arbitrary radial function that defines a set of gauge transformations. Note that 
1^ = gives back the original regular Hamiltonian, while ^ = ^1 together with = /c = 1 at 
large distances yields the return string configuration considered in ref. [271]. Thus the interpolation 
between regular and singular behavior is accomplished by the boundary conditions ^(0) = and 
^(00) = .^1. The transformed Dirac Hamiltonian becomes 



H 





crp 



^"^p 



cr • ^ 



p \cr ■ if 



5. 



int : 



cos(A) 



1 

-1 



+ i sin(A) 



1 

-1 



n • T 



+ 



2dp 



-a- ■ p 



(T ■ p 
-a ■ p 



n • T 



n / —a ■ (p (7 ■ ip 
2p \ cr ■ (f -a ■ (f 



/Gsin(A)/G(A) + ifc - l)sin(0/G(-e) 



(10) 



The new gauge function ^{p) is hidden in the difference A(p) = ^1 — (,{p) which appears both 
explicitly and as the argument of the space-dependent weak isospin matrix 



Icix) 



' — sin(x) — zcos(2;)e 
i cos(a;) e~*"'^ sin(3;) 



(11) 



All explicit matrices in eq. (jlOp act in spinor space. Together with the boundary conditions for the 
string profiles /c and eq. (jlOp defines a well-behaved scattering problem for which a scattering 
matrix and, more generally, a Jost function can be straightforwardly computed. Moreover, the 
Born series to these scattering data can be constructed simply by iterating -ffint- 

We will renormalize the calculation by subtracting orders of the Born series and adding these 
contributions back as the corresponding Feynman diagrams. It should be mentioned that although 
the Jost function is gauge invariant, neither the Born series nor the individual Feynman diagrams 
associated with eq. (jlOp is gauge invariant, and so the Born subtracted phase shifts or Jost func- 
tions will also depend on the gauge. That is, these quantities are functionals of ^(p). However, 
the gauge-dependent terms subtracted from the phase shifts correspond exactly to the gauge- 
dependent finite parts in the Feynman diagrams, while the counterterms, which parameterize the 
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ultraviolet singularities, are gauge-independent. The net effect is that individual pieces of the 
spectral approach to the vacuum energy will be gauge-dependent, but the combined expression is 
not. 

The formulation in eq. (|10p describes the physical string without the need to introduce artificial 
return strings to unwind the topology. In particular, the (tedious) separation of the return string 
contribution from the bound state spectrum of the physical string is no longer required. Moreover, 
the gauge function ^{p) can be taken to have support at a moderate distance, so that there is no 
need for non-trivial fields at very large radii and, as a consequence, the need for extremely large 
angular momenta is avoided. 

To solve the scattering problem for the Hamiltonian, eq. (|10p . in two space dimensions, we first 
introduce grand-spin states to take care of the angular dependence. For a fixed angular momentum 
i, there are four grand-spin states, characterized by the quantum numbers ±1/2 for spin S and 
isospin /, 



{^■,SI\i + +) =e 



i{£+n)ip 



{^■,SI\e-+) = ie*(^+"+i)^ 



{v;SI\i + 



{v;si\e 



-I e 



(12) 



The angular dependence is thus separated from the radial dependence by the ansatz 

^iip,ip)= Yl {{p\{^; si\)\eesj). 

For each value of the angular momentum I, this decomposition turns the Dirac equation 



(13) 



(14) 



with the Hamiltonian given in eq. (jlOp . into a 8 x 8 system of ordinary first order differential 
equation for the radial functions in the spinor states 



{p\ei + +) 



{P\ee-+) 



9i{p)\i-+) 

f3{p)\e-+) 
93{p)\i + +) 



{p\ei + 



f2ip)\e + 
92{p)\i- 

h{p)\e- 
9Aip)\i + 



(15) 



where we have suppressed the energy label (e) on the radial functions. It is convenient to combine 
these eight functions in a vector notation 



/ 



/2(P) 

fsip) 

\h{p)J 



and 



(9i{p)\ 
92{p) 
93{p) 

\94{p)J 



(16) 



In terms of these vectors, the static Dirac equation in each angular momentum channel takes the 
form of two coupled real 4x4 systems, 

(e -m)f = Vuu / + {-CDu + Kd) 9 



{e + m)g = {CDd + Vdu) / + Vm 9- (17) 
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The 4x4 matrix C = diag(— 1, — 1, +1, +1) is constant while and contain the radial 
derivative operator as well as the angular barrier terms. The coupling to the background profiles 
of the boson fields emerges via the matrices Vij. Detailed expressions for D^, and Vij are 
listed in appendix |Bl The ODE system eq. (flTl) is the basis of the spectral approach to the string 
problem. 

For the gauge profile £,{p), any smooth function with ^(0) = and ^(oo) = will do. For 
simplicity and to avoid possible singularities at /> — )■ 0, we choose again a Gaufiian profile 

= 6 [l - exp i-p'/wl) ] (18) 

with a new width parameter w^. As explained earlier, the scattering matrix without Born sub- 
tractions and the complete quantum energy should be independent of the choice of gauge and 
thus inde pen dent of the width parameter w^. This has been verified numerically to a fairly high 
precision [26(. 



IV. SPECTRAL METHOD 

In this section, we present the details of our approach to compute the fermion contribution to the 
vacuum energy of the string. To make the exposition clearer, we have moved overly complicated 
expressions and all technical derivations to the appendices. However, the complete method is still 
quite involved due to the many contributions that enter. We will continue the discussion of the 
variational approach for charged strings in section |V] and present numerical results in section IVll 
The calculation of the fermion quantum energy is based on the Dirac equation (I14p . From the 
solutions to this equation we infer a number of distinct contributions to the energy of the string, 

Ef = Es + ^FD + (19) 

In physical terms, these three contributions are 

Es: the non-perturbative vacuum polarization due to the string background, with the divergent 
low-order Feynman diagrams taken out by subtracting leading terms in the Born expansion. 
This piece also includes the bound state contribution to the fermion determinant; 

E-pj)-. the perturbative contribution of the low-order Feynman diagrams to the vacuum polarization 
energy, combined with the counterterms for proper renormalization. This compensates for 
the part that has been taken out of Es by means of the corresponding Born expansion; 

E]^: the binding energy due to the single particle bound states that are explicitly occupied to give 
the string a fermion charge Q. More precisely, E\^ = X] ~ Q™- measures the energy 

occ bs 

of the populated levels relative to the same number of free fermions. We will describe this 
contribution in the next section. 

Each of these pieces is separately finite; the first two terms are not gauge invariant, but their sum 
is, and so is E\^. 

In this section we focus on the renormalized vacuum polarization energy 

Eq = Es + Efy) . (20) 



Potential ambiguities in Eg that could originate from the ultra-violet divergences are fully removed 
by the identification of terms in the Born series with Feynman diagrams. The most important 
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feature of Eg is the possibility to impose renormalization conditions from the perturbative sector 
(MS or onshell), although the calculation is completely non-perturbative, including all orders 
in i?int- We then combine Eg with the classical energy E^y required to form the bosonic string 
background. Quite generally, we expect E^i + Eg > once quantum fluctuations are included, 
since otherwise we would have an instability of the true vacuum to cosmic string condensation, 
which should obviously not happen. 

In the following subsections, we will give brief accounts for each contribution to eq. ()20|) . More 
details can be found in the appendices. 



A. Jost function and Born subtractions 

The Born-subtracted vacuum polarization energy Eg has contributions from bound and scat- 
tering states. These two contributions are combined in the Jost function for imaginary mo- 



menta [29l . |34| |. To compute Eg it is therefore sufficient to solve the scattering problem as in 



ref. 23|: For every energy |e| > m, the fermion system eq. (I17p has eight real linear independent 



solutions {f,g)- In the case without a string background, these solutions are Bessel functions 
of integer order with the argument z = kp, where k = \J — ra^ > 0. Instead of taking the 
(real) regular and singular Bessel functions Jjy(z) and Yy{z), respectively, we can formally let (/,5) 
have complex coefficients and take Hankel function solutions instead. In this case, both the real 
and imaginary parts of {f-,g) are (linearly independent) solutions, or equivalently {f,g) and their 
complex conjugates are independent solutions. 

To describe the coupled channel scattering problem, it is convenient to put the four free linearly 
independent complex solutions for / and g onto the diagonal of two 4x4 matrices 

= diag(//];)„(M,^f (M,^£l+i(M,^il\(M) (21) 

Ud = diag(F«„^,(M,^ll\(M,^l+l(M,^f (M) , (22) 
which describe out-going asymptotic fields since 

F(i)(z) = Uz) + iY,{z) ,/Ae<--f--|-) , (23) 

V vrz 

as z — )• oo. With this notation, the j*^ linear independent solution is 

where denotes the j*^ row of the matrix T-L. For convenience we omit the orbital angular 
momentum index By construction, the complex conjugate matrices, 7^* ^ describe incoming 
spherical waves. Furthermore, we have defined the relative weight of upper and lower Dirac com- 
ponents as 

^ ^ - (24) 



e + m k 

For later analytic continuation we must ensure that the phase of the Jost function is odd for real 
momenta under k — t- — fc, which requires the branch cut structure of the square root be defined 
using either of the two expressions listed above. 

To describe the coupling of the four channels in the actual scattering problem, it is convenient 
to put the linearly independent solutions for / and g again in the rows of a 4 x 4 matrix, and factor 
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out the free part to get simple Jost boundary conditions, 

(/), [J-.H^. and 

(25) 

ig)j K [g-n,]^ and k [G* ■ K]^ . 

We substitute these ansdtze into the Dirac equation, eq. (fT7|) . and find first order differential 
equations for the matrices J- and Q. This is explicitly carried out in appendix IB 1[ The solutions 
to eq. ()B6P with the Jost boundary conditions 

lim T{p) = lim g{p) = l (26) 

p— >oo p— >oo 



define Jost solutions to the initial Dirac problem via the representation in eq. (j25p . The physical 
scattering solution is the linear combination which at large distances is the superposition of incom- 
ing and outgoing free spherical waves and obeys the regularity condition at the origin. The relative 
weight of the incoming and outgoing waves defines the scattering matrix S. Hence the physical 
scattering solution for the J^-type (upper) components reads 

^> = T* ■ni + {T-nn)-S. (27) 

The corresponding ^-type (lower) components are obtained by replacing T ^ Q and T-Lu ^T-Ld- 
The physical scattering solution must be regular at the origin p = 0. From this condition we 
extract the scattering matrix in one of two equivalent ways 

S = -limn;:'-T~'-T*-n: (28) 
S = - lim • • • . (29) 

The phase convention in eq. (|27p is chosen to reproduce 5 = 1 for the non-interacting case which 
has J^^^^ = g^'^^ = 1. The equality of the two representations from the system of coupled differential 
equations is a good check on our numerics, as is the requirement that S be unitary. 

It should finally be noted that all the interaction matrices eq. (IB9[) are linear in the background 
profiles eq. (|B4p . so that the ODE system for the Born approximation can simply be obtained by 
iteration with the ansatz 

oo oo 
-^Born(p) = 5^-^(*Hp) and gBorn(p) = ^a(*)(/5), (30) 

i=0 i=0 

where the superscript denotes the order of the interaction Hamiltonian Hi^^t in eq. (jlUp . The zeroth 
order solutions are J^^^\p) = Q^'^^p) = 1 and all subsequent contributions are subject to the 
boundary conditions 

lim -F(*)(p) = lim g'^^{p) = 0, i = 1,2,3... 

p— >oo p— >oo 

The explicit form of the iterated system of differential equations for the Born approximations of 
order i = 1 and i = 2 can be found in appendix |Bj Though the i = 3 and i = 4 orders also yield 
divergences, we do not discuss them explicitly because we employ a numerically less costly method 
to handle these logarithmic divergences, as described below. 



11 



B. Interface formalism 



The 4x4 scattering matrix S derived in the last subsection yields the four eigenphase shifts and 
thus the shift in the two-dimensional density of states j35l | 



d6e 



dk 



i d 
2tt dk 



lndet5£(A;) , 



(31) 



where the sum runs over the four scattering channels for a given grand spin channel, which we 
label by the associated orbital angular momentum i. To turn this two-dimensional density into a 
three-dimensional energy (or energy per unit length of the string), we have to deal with the trivial 
dynamics along the string symmetry axis of the string. This is a typical application of the interface 
formalism developed in ref. [s^]- The modifications of the usual spectral method are simple: 

1. The integration over the momentum conjugate to the coordinate of translational invariance 
remains finite due to sum rules for scattering data 3l|, [S^j that are generalizations of Levin- 
son's theorem. 



2. When integrating over momentum k, the density, eq. (I3ip . must be multiplied by a kinematic 
factor that differs from the usual one-particle energy e = \/'m'^ + k'^. 

3. More Born subtractions are required to make the momentum integral and angular momentum 
sum convergent. This corresponds to the larger number of divergent Feynman diagrams in 
three dimensions. 



Then the interface formula for the vacuum polarization energy per unit length of the string is 



E 



(N) 



iE °' 



°^ dk 





+ ?7T,^)ln 



k^ + 



dk 



N 



+E 



(gj- ^) (ej- ^) + 



(32) 



where the notation [•••Jat refers to the quantity in the brackets with its first N terms of the 
Born series subtracted. For the string problem in three space dimensions, we need > 4 to 
ensure convergence of the momentum integral although, as described below, we will use a different 
subtraction in place of the = 3 and = 4 cases. Here, e^^^ gives the energy of the j^^ bound 
state in the angular momentum channel I, and is the degeneracy in that channel. For the string 
background, we have 



Do 



> 



-n 
-n . 



(33) 



where n = 1 is the Higgs winding number introduced in the string configuration of eqs. ([5]) and ([6]). 
The renormalization scale /i emerged from the integration under item 1. It cancels due to the 
same sum rules. For convenience we usually set ^ = m. The phase shifts can be extracted from 
the scattering matrix or, equivalently, from the Jost-like matrices J- and Q introduced in the last 
subsection, 



5Ak) = - Indet lim Fi[p, k) ^TUp, k) = - Indet lim Qi{p, k) ^Gpip, k) 

I p^O I p^O 



(34) 
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where we have restored all the arguments. In deriving eq. (j34p from eq. (|29p we have used the 
cyclic property of the trace and the fact that as p — t- the Hankel functions are dominated by their 
imaginary parts. 

As indicated in the previous subsection, it is convenient to evaluate the expression (j32p in the 
complex /c-plane because after rotating to the imaginary axis the explicit bound state contribution 
is automatically canceled by the pole contribution from Cauchy's theorem, leaving only a single 
integral along the cut on the positive imaginary axis 2^, 34 1. There is another important technical 



reason to rotate to imaginary momentum. We need to sum over angular momentum i and inte- 
grate over radial momentum k after subtracting sufficiently many terms of the Born series. This 
procedure is numerically cumbersome because these functions oscillate in k, which can make it 



impossible to exchange the sum and integral 36|] because they are not absolutely convergent. This 



obstacle is also avoided by analytically continuing to imaginar y m omenta t = ik and performing 



the integrals in the complex plane along the branch cut t > m j36l |. 

The analytic continuation for the Dirac equation is conceptually different from the well-studied 
Schrodinger case because e = ±Vk'^ + m'^ causes the complex momentum plane to have two sheets. 
So on the real axis we have to pick one sign, continue to complex momenta and compute the Jost 
function on the imaginary axis. This procedure must then be repeated for the other sign and 
then all discontinuities must be collected at the end. In the present problem we are fortunate 
because the solutions to the Dirac equation exhibit charge conjugation symmetry along the real 
axis. Therefore det(5) does not change under e — )• — e and there is no additional discontinuity 
in the Jost function choosing either sign. Moreover, the Jost function is real on the imaginary 
axis, as in the Schrodinger problem. However, the way this comes about in the string problem 
requires us to be careful when constructing the Jost function for complex momenta. This procedure 
is described in appendix |B] and results in the replacement of the phase shift 6{k) (and its Born 
expansion) by i'(t), the (modified) logarithmic Jost function for imaginary momentum. For this 
to work it is essential to have k odd under sign reflection of real k. The resulting Jost function 
itself is a continuous function in the upper complex momentum plane and the branch cuts in the 
Dirac equation do not carry over to i'{t). The only discontinuity arises from the logarithm under 
the integral in eq. ()32p . which is 27r. Finally an integration by parts yields a simple expression for 
the Born subtracted vacuum polarization energy. 



-—jdttY,D,Mt)]N- (35) 



Here we have interchanged the integral with the angular momentum sum, which is possible on the 
imaginary axis [s^ . After a final change of variables t — t- r = Vf^ — rn?, we obtain eventually 

= / E MVr^+^^) ■ (36) 

Eq. (j36p is our master formula for the phase shift contribution to the vacuum polarization energy 
per unit length of the string. 

C. Feynman diagrams 

The Born subtractions in the integrand of eq. (|36p must be added back in as Feynman diagrams. 
The latter are most easily derived by expanding the fermion determinant representation of the 
(unrenormalized) vacuum polarization energy. 



A = -TLEq = -i In det 



m + Hi 



(37) 
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Both the time interval T and the length L of the string factorize as T — )■ oo and L — )• oo because 
the string background is static and translationally invariant. The interaction part of the Dirac 
operator can be separated in various spin structures, 



Hi = L^-f''PL + h + ipj5 
where the fields L^, h and p are isospin 2x2 matrices, 



(38) 



Lo = 0, 

h{p) = -anip) 1 
The isospin matrices Iq, defined in eq. (jlip . and 

Ip{ip) = n • T = 



Up) = 2 ar{p) Ip{^) p + 2 [acip) /g(6 - ^P)) + Ici-ap)) 
p{p) = -ap{p)Ip{(p) . 



e*"'^ 



(39) 
(40) 



contain the entire dependence on the azimuthal angle (p. The profile functions fnip), fdp) and 
^(p), cf. eqs. (I8|) and (jlSp . determine the radial behavior of the coefficient functions ar{p), Oip{p) 
and a^{p). They are explicitly listed in eq. ()B4p of the appendix. The Feynman series for the 
effective fermion action (determinant) is now 



A = -TLEq = -i In det {i^- m) + ^ ^ Ij Tr {ip 



N 



N=l 



N 



m) 



Hi 



N 



(41) 



where the first term corresponds to the free vacuum energy without a string background. It is 
automatically removed in the spectral method by the difference in eq. (j3ip . For fermions in three 
dimension, all diagrams up through A'^ = 4 are divergent and thus subject to renormalization. 
While the calculation of the corresponding Born subtractions up to fourth order is not particularly 
hard, the evaluation of the higher order Feynman diagrams with up to four nested Feynman 
parameter integrals and an equal number of Fourier transformations of the string background is 
very cumbersome. A better approach is the so-called fake boson method introduced in ref. 
which we will describe next. 



37|, 



D. Fake boson approach and renormalization 

The first and second order fermion Feynman diagrams contain both quadratic and subleading 
linear and logarithmic ultra-violet divergences, so that a precise identification of the terms in the 
Born expansion with the Feynman diagrams must be made separately for each term in the angular 
momentum sum. On the other hand, the third and fourth order fermion Feynman diagrams only 
cause logarithmic divergences, which are much easier to cope with, because the sum 

5(t) = 5;Z),h(t)]2 (42) 
e 

is finite. However, after multiplication by t, the integral in eq. (I35p is logarithmicly divergent. So 
instead of subtracting the complete third and fourth order terms in Hi from the sum in eq. (j42p . it 
is sufficient to just subtract any function As(t) of momentum with the same ultra-violet behavior, 
provided that the following conditions are met: 

1. the subtraction As{t) should have the same analytic properties with respect to complex 
momentum arguments; 
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2. formally its contribution to the vacuum polarization should be identifiable as a Feynman 
diagram that can be combined with the available counterterms 



C^t = ci tr(G^^ G^^) + C2 tr [ [D^^]^ D^^'\ + cg [tr(«>t$) - 2v'^'\ + C4 [tr($t$) _ 2v 



(43) 



to cancel all ultra-violet divergences. The perfect candidate is the second order contribution from 
a boson scattering off a radially symmetric potential V{p)- Prom the properties of the (bosonic) 
scattering problem [sH], we know that its Jost function has the required analytical properties and 
its contribution to the vacuum polarization energy can be expressed as a (very simple) Feynman 
diagram. It only remains to adjust its strength to accomplish the required subtraction. This 
fake boson scattering problem also has a partial wave decomposition and we subtract the sum of 
the logarithm of the second order fake boson Jost function from the sum in eq. (j42p . Since the 
subtraction is not carried out channel by channel, the exchange of i-sum and i-integral is crucial 
for this approach to work. 

To describe the method in detail we define An to be the contribution of order (-ff/)" in the sum 
in eq. ([iT]) . 

1. The first order diagram = 1 is linear in the interaction Hj and local, including all fi- 
nite parts. Thus the entire diagram is proportional to the spacetime integral of the C3- 
counterterm in eq. (j43p . We fix the corresponding counterterm by the no-tadpole condition 

Ai = 0, which ensures that the vev of the Higgs field is kept at its classical value v. This 
condition completely fixes both the divergence and the finite part in the cs-counterterm, cf. 
eq. (ICT6]) . 

2. The second order diagrams N = 2 give contributions to the various propagators whose con- 
tributions to the vacuum polarization energy are quadratically divergent at large momenta, 
for which a careful regularization is required. Due to gauge invariance, the coefficients ci, C2 
and C4 in the counterterms in eq. ()43p can unambiguously be determined by the two-point 
functions that emerge at order N = 2. Hence we do not need to compute the full Feynman 
diagrams at orders = 3 and = 4. 

3. Although we do not need the full diagrams, we do need to precisely subtract the divergences 
from ^3 and A4. In dimensional regularization (D — )• 4), these logarithmicly divergent pieces 
read 

4-D r d^k 



AOO 

-^3,4 



(2vr) 



D 



[k^ - 1 + ie)' 



(44) 



where T and L are the (infinite) lengths of the time and z-axis intervals, respectively, and cp 
is a complicated integral over the radial profile functions, cf. eq. (jClOp . The key observation 
for the implementation of the fake boson approach is that the divergence in eq. (j44p is also 
contained in the two-point function of a simple scalar field that fluctuates in a (flctitious) 
background potential V{p)- In fact, the divergence in the second order boson diagram has 
the form of eq. (j44p with cp replaced by 



1 



CB 



dppV{pf. 



(45) 



By properly scaling V{p) with ^JcfJcb., we can match the divergences from eq. (I44p . The 
equivalence of the Feynman and Born expansions implies that the combination of s(t) with 



As(t) 



-(2 

CB ^ 



it) 



(46) 
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is finite when integrated according to equation (f35l) . Here (t) is the second order Born 
approximation for logarithm of the Jost function on the imaginary axis in the fake boson 
problem and the associated degeneracy factor in the partial wave decomposition is Di = 
2 - 5o,£- 

4. The subtraction in eq. ()46p must be compensated by adding the corresponding second order 
fake boson diagram. Since the divergences of the fake boson and the fermion problem have 
been carefully matched, the fermion counterterms from eq. (j43p are sufficient to render the 
relevant fake boson diagram finite. As a consequence, only the renormalized fake boson 
diagram must be added back in, cf. appendix IC 2[ We are then fully prepared to compute 
the vacuum polarization energy in any renormalization scheme. We first consider the MS 
scheme, which is defined by setting = in the counterterm coefficients 




In this scheme the dependence on the model parameters is simple. The computational 
advantages of first considering the MS scheme will be discussed thoroughly in section IVII 

Let us summarize the result in the MS scheme and carefully describe the angular momentum sums. 
First we construct the subtracted logarithmic Jost function for imaginary momenta 

^max -^max 

Kt)= lim V I)£[z^Ki)]2 + — lim Y.^iVf\t) . (48) 

From eq. (08]), we can compute the phase shift contribution to the vacuum polarization energy, 

Es = / dTri/(Vr2 + m2) . (49) 

27r Jo 

The complete vacuum polarization energy in the MS scheme is then the sum 

E^ = Es + AEfd , (50) 

where AE'fd = ^-^fd ~^ ^-^B is the sum of the renormalized values (finite parts in MS) of the 
second order fermion and fake boson diagram. Explicit expressions for these contributions can be 
found in eqs. ()C8p and (jClSh . As a further test of the approach we verify numerically that E^j^ 
remains unchanged when the boson potential V{p) is modified. 

To make contact with the electroweak theory, it is convenient to re-adjust the finite pieces in 
the counterterms such that they match the so-called on-shell scheme. In addition to the already 
implemented no-tadpole condition that fixes C3, we thus require 

• The pole of the Higgs propagator remains at the tree level mass, m/j = m^^^ , with unit residue. 
This fixes the coefficients C2 and C4 and ensures the usual one-particle interpretation of the 
states created by the asymptotic Higgs field. 

• The residue of the gauge field propagator (in unitary gauge) is unity, so that asymptotic 
VF-fields create one-particle -boson states. This condition determines ci. 

The position of the pole in the gauge boson propagator is then a prediction, i.e. the physical 
VF-boson mass receives radiative corrections. These corrections are determined by the implicit 
solution to eq. ()C17p presented in appendix IC 31 
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The on-shell and MS schemes are related by finite changes of the counterterm coefficients, 
so that Cj 7^ in the on-shell scheme; explicit expressions can again be found in appendix IC 31 
The modification to the vacuum polarization energy eq. ()50p due to the change in renormalization 
scheme is then simply the energy from the counterterm Lagrangian eq. (|43p . with the coefficients 
Cj replaced by their finite pieces q listed in eq. (jClOp . Since the counterterms are local, this modi- 
fication amounts to a radial integral similar to the classical energy in eq. ([7]), which is numerically 
inexpensive. Hence the vacuum polarization energy in the on-shell scheme is 



En 



+ EcT = Es + A^D + ^^B + EcT 



(51) 



The explicit expression for the counterterm contribution reads 



EcT 



27r 



J pdpjsin^^i 



n 



C2V^fH{l-fGy 



4C1 

9' ^' 



(52) 



The counterterm coefficient C3 does not appear explicitly because this counterterm receives no 
correction in passing between the MS and on-shell schemes. 



V. CHARGED STRING 

As already discussed in refs. 2^, 27 1, the fermion vacuum energy is negative for narrow strings 
and thus provides some binding. However, for physically relevant model parameters, eq. ([4]), 
it is insufficient to overcome the large classical energy. The central mechanism for overcoming 
the classical energy cost is to populate the numerous fermion bound states that emerge in the 
background of the string, which as a result assigns charge Q to the string. If the energy of 
equally many free fermions Qm is larger than the total energy of the string (the classical, vacuum 
polarization and contribution from populated levels combined), we have succeeded in constructing 
a stable charged string. Quantitatively, this requirement corresponds to E^i + Ef < 0, cf. eqs. d?]) 
and ()19p . It prevents the direct decay into fermions and only leaves charge non-conserving decay 
channels, where the decay rate is heavily suppressed due to the sphaleron barrier. The direct decay 
into lighter fermion doublets is also suppressed, since we do not have fiavor mixing in our model. 
To carry out this procedure, we first need to find the bound state energies in the string background. 



A. Bound states and box diagonalization 

Since our system is translation invariant in the z direction, we begin by finding Dirac bound states, 
€i, of the two-dimensional problem, eq. (jl4p . Each bound state we find will then correspond to a 
family of bound states in the three-dimensional problem, indexed by the transverse momentum ■ 
We carry out the two-dimensional bound state calculation by putting the string in a large cylin- 
drical box of radius R ^> and imposing the boundary condition that no net fiux runs through 
the surface of the cylinder. This boundary condition discretizes the possible radial momenta in 
each angular momentum channel through the roots of certain Bessel functions, cf. eq. ()A2p for the 
case of unit winding, n = 1. We can thus take a countable set of grand spin solutions, eq. ()12p . to 
the free Dirac equation and express the fully interacting string Hamiltonian as an infinite matrix 
in this basis. The relevant matrix elements are again presented in appendix [Aj Upon truncating 
the set of free solutions by including an effective UV cutoff A on the discrete momenta, we are 
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thus left with a large matrix diagonalization in each grand spin channel to determine the fermion 
eigenstates in the string background. Typical matrix sizes are (1600 x 1600) including Dirac indices 
We then find the energy eigenvalues numerically by diagonalization. In the finite box, of course, 
all energy levels are discrete and there are no continuum states. In the limit ii — )■ oo and A — )■ oo, 
the highest energy levels in the quasi-continuum will still fluctuate considerably, but the low-lying 
bound state spectrum of states with energy smaller than m, which become bound states in the 
i? — 7- oo limit, should remain stable. This was indeed observed for moderate values, A ~ 8m and 
R ^ 75/m. 

It should be noted that bound states occur predominantly in the lower angular momentum 
channels, as we would expect since the higher channels contain an increasingly large centrifugal 
barrier. Depending on the width of the background profile, we see bound states in as many as 10 
channels, or only in the single channel i = —n = —1, which is the channel that contains an exact 
zero mode for = '^/2. 

A good numerical test on our diagonalization procedure is the gauge invariance of the Dirac 
Hamiltonian and thus of the bound state spectrum. In our specific case, this means that the 
low-lying bound state energies must remain constant when the gauge transformation profile ^(p) 
is modified. We have confirmed this behavior for simple scale and width changes in ^(p). 



B. Populating the Bound States 

Having determined the set of bound state energies in the two-dimensional problem, we now 
integrate these results into the full three-dimensional calculation. Let < < m represent the 
energy of one of these two-dimensional bound states. In the full three-dimensional problem, we 
will then have a family of bound states with energies [e? H-p^] ■ 

For a given charge Q, each of these families of bound states will be filled up to a common 
chemical potential^ A^(Q) ^ '"^j to minimize their contribution to the energy. If the towers of states 
built upon two different had different upper limits, the energy would be lowered by moving a 
state from the tower with the larger limit to that with the lower one, without changing the charge. 
Since states with [ef +p'i]^^'^ < n are filled while states with [ef > remain empty, we have 

a Fermi momentum Pi{^x) = [/x^ — ef]^''^ for each bound state. By the Pauli exclusion principle we 
can occupy each state only once, and so we find the charge density per unit length of the string 

Q(p) = i^P,(/i), (53) 

where the sum runs over all bound states available for a given chemical potential,^ < p. Of 
course, this sum involves different partial waves, so we have to include the corresponding degeneracy 
factors. 

Eq. (j53p can be inverted to give p = fJ,{Q). In numerical computations we prescribe the left- 
hand-side of eq. (j53|) and increase /i from min{|ei|} until the right -hand-side matches. From this 
value fi = fi{Q), the binding energy per unit length 



1 f^^W 

'b(Q) = - 2^ / dp^ 

^ Jo 



^ In what follows the chemical potential fi should not be confused with the redundant scale introduced in eq. (|32p . 
^ Ambiguities in this relation due to different boundary conditions at the end of the string show up at subleading 
order in 1/L, where L is the length of the string, and can thus be safely ignored. 
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ei<fi 



P,(/x)(//-2m) + e2 In 



(54) 



can be computed as a function of the prescribed charge. In this manner the total energy becomes 
a function of the charge density of the string. In our search for a stable string, then, we specify 
the charge Q, and, among background configurations with sufficient binding to accommodate this 
charge, we vary the ansatz parameters to minimize the total energy to see if we find a bound 
configuration. 



VI. NUMERICAL RESULTS 

In this section we present the numerical results combining all the contributions to the string 
energy in our variational ansatz. We measure the variational parameters wh and wq in inverse 
fermion masses. The dimensionless vacuum polarization energy per unit length Eq/vn? then does 
not explicitly depend on the coupling constants / and g in the MS scheme, and depends only 
weakly on these constants in the physical on-shell scheme through the logarithmic dependence 
introduced by the renormalization conditions. This property simplifies the numerical analysis 
because then their variation solely affects the classical and counterterm energies, both of which are 
local functionals of the profile functions and hence easy to compute. 



We have already presented first results for the vacuum polarization energy, eq. (jSip in ref. [26|. 
In particular, we have verified our results numerically by checking that they are independent of the 
shape of the gauge function ^ (p) . This result is a consequence of gauge invariance, but it is nontrivial 
because the individual Born terms and Feynman diagrams are not explicitly gauge invariant — 
only the combination of all of them is. As a result, this invariance verifies the equivalence between 
the Feynman diagram contribution and the Born subtractions (including the fake boson part) in 
eq. ()5ip , which is central to the application of spectral methods in quantum field theory [29[] . 

The computation of is numerically most costly. The main reason is that we have to go to very 
high angular momenta in the sum in eq. (f48]) . Typical values are £max = 500, . . . , 800 depending 
on the width of the background field. To capture the behavior of the integrand in eq. (j49p we 
consider about 40 points in the interval < r < 8. Since the integrand of Es does not oscillate 
when computed from imaginary momenta, we can accurately estimate the contribution from r > 8 
from an inverse power-law behavior. 

In figure [2] we show the result of this numerical computation. The wider the background 
fields, the weaker the dependence on the angle that parameterizes the gauge boson contribution. 
Surprisingly, we see that the vacuum polarization per unit length is quite small. Even for large 
widths it does not exceed a fraction of the fermion mass squared. With the exception of very small 
widths, the vacuum polarization turns out to be positive. Hence there is no indication that the 
vacuum polarization energy from the fermions can stabilize cosmic strings since the classical energy 
is larger by orders of magnitude, unless the coupling constants are /, 5 ~ 0(10). For example, 
see figureO which shows the classical energy for the standard model parameters, which are 0{1). 
The derivative terms of the classical energy decrease quadratically with / and g while the Higgs 
potential decreases like 1//^ for fixed Higgs mass. As a result, increasing the coupling constants 
could lead to binding for thin strings, but such configurations contain large Fourier components, 
which for /, 5 ~ 0(10) reach the vicinity of the Landau ghost pole. Hence any such binding is 
obscured by the existence of the Landau ghost, c/. appendix C.4, which arises when including 
quantum corrections in a manner that does not reflect asymptotic freedom. Here it is due to 
the omission of quantum corrections from fluctuating gauge boson fields. The estimate for the 
Landau ghost contribution discussed in the appendix suggests that the issue can be safely ignored 
for f,g<5. 
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FIG. 2: Vacuum polarization energy as function of the angle for different values of the width parameters 
wh and wq in the on-shell renormalization scheme. The physically motivated model parameters, eq. ([4]), are 
used. The dots refer to actual computations, while the lines stem from a cubic spline. We also show the results 
obtained for the fit to the Nielsen-Olesen profiles, c/. figure [1] These results do not include the combinatoric 
color factor Nc- 
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FIG. 3: Classical energy for the standard model parameters as function of the ansatz parameters wh, wg and 
^1, including the fit to the Nielsen-Olesen profiles. The model parameters are again from eq. ([4]). 



Gradient expansions j38l | for quantum mechanical expectation values suggest that the energy 
gain from populating bound states can be estimated from a spatial integral over some fractional 
power of the potential in the wave-equation. Scaling arguments show that the energy from the 
populated bound states increases quadratically with the width parameters wh and/or wq, regard- 
less of the specific power in the expansion^. Since also the dominating classical energy increases 
quadratically with wh (from the Higgs potential), populating the bound states might balance the 
large classical energy already at small coupling constants and moderate widths, because the Higgs 
potential scales like 1//^ when the Higgs mass is fixed while the fermion contribution is not sen- 
sitive to any change in / when we scale our definitions of physical quantities with the fermion 
mass m = vf. Hence our strategy to construct a stable string type configuration is to balance 
the classical boson energy with the fermion quantum correction by considering wide strings and 
increasing the Yukawa coupling. Simultaneously we must keep fixed the charge associated with 
populating the fermion bound states. 

Before we will consider the total energy we would like to discuss the fermion part, Eq + E^. In 



^ More precisely, the energy gain involves both the summed energy eigenvalues and the charge. Both can be expressed 
by such integrals with different powers, though. 
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figure m we show the total fermion energy of eq. (I19p . as a function of the charge density per unit 
length of the string. As described in the previous section, we can compute the binding energy as a 
function of the charge for a given background configuration. After adding the vacuum polarization 
energy computed for that background, we get the parabolic curves in figureHl These lines terminate 
at the point where all available bound states are populated. We then search for the configuration 
that minimizes the energy. For small charges, we obtain thin strings, while larger charges lead to 
wide strings, as shown in figure |4l Surprisingly, the resulting envelope that describes the minimal 
fermion energy as a function of the charge density is a straight line with (approximately) vanishing 
y-intercept. This straight line stems from a delicate balance between the vacuum polarization 
and binding energies. Because this extrapolation yields a vanishing y-intercept, we deduce that 
very narrow strings have vanishing vacuum polarization energy. This interpolation overcomes the 
Landau ghost problem of the direct calculation. 

From several hundred configurations for which we have computed both the vacuum polarization 
energy and constructed the bound states, we identify the one that minimizes the total binding 
energy 

^tot = E^i + E^ + Ek (55) 

for a prescribed charge. The corresponding result for the minimal binding energy is displayed in 
figure El We can see that the optimal binding grows linearly with Q. The steep slope at very 
small charges is an artifact of restricting our ansatz to configurations with wh,wg > 2 to avoid 
unphysical effects from the Landau ghost. 

As mentioned above, we increase the Yukawa coupling from its top-quark motivated value 
/ = 0.99, while all other model parameters are taken from eq. ([4]). Increasing the charge also 
increases the width of the optimal string. For / ~ 1.6 the classical and quantum contributions 
balance and the total energy is essentially independent of the width. Increasing the coupling further 
yields a negative energy (in comparison to equally many free fermions) and stable configurations 
exist. Not surprisingly, the minimal charge for which there are stable configurations decreases 
quickly as / increases. For / = 1.7 it is Qmm ~ 10m = 17v, while for / = 1.9 stable configurations 
exist already at Qmin ~ 3m = 5.7 v. 

We next discuss the structure of the stabilizing configuration. We find that the fermionic part of 
the binding energy is insensitive to the angle ^i, as shown in figure [6j As a result, the dependence 
of the total binding energy on stems entirely from the classical part, which is clearly minimized 
for ^1 ~ 0, since in that case the the gauge fields vanish and we have only a charged Higgs field. 
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FIG. 5: (Color online) Total energy as a function of the charge for various values of the Yukawa coupling constant. 
The color degeneracy Nc = 3 is included. 
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FIG. 6: (Color online) Fermionic contribution to the string binding energy per unit length as a function of charge 
density per unit length, in units of the fermion mass, for a variety of values of with wh — 6.0 and wg ~ 6.0. 

with only the non-diagonal elements in eq. ^ differing from zero. 

In figure [7] we display the chemical potential that minimizes the binding energy for a prescribed 
charge. Its construction is discussed in section|Vl The cusps arise because our sample configurations 
are not continuous in the variational parameters. As we increase the charge, the minimizing 
configuration jumps among these possibilities. 

The strong deviation from fx = m at low charges (where the various graphs overlap) is again an 
artifact of not considering very narrow string configurations. We see that at the limit of binding 
(/ = 1.6) almost all bound states are populated. As the binding increases, the chemical potential 
decreases, leaving the states below threshold un-occupied. 

We have seen above that binding increases with the Yukawa coupling. Table |T] indicates that at 
the same time the profile functions get wider, while the critical charge at which binding sets in (i.e. 
£'tot(Qmin) = 0) decreases with the Yukawa coupling. As a result, the width of the critical profile 
actually decreases. We find these widths to be 5.5/m = 3.2/v, 4.0/m = 2.2/v and 3.5/m = 1.8/v 
for / = 1.7, 1.8 and 1.9, respectively, and the typical extension of a bound charged string is about 
0.003fm. 

Finally let us estimate the total mass of the bound string. In the regime where it is only 
slightly bound, we have Stot ^ Qm. Typically we observe binding for Q ~ 5m. Hence a reasonable 
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FIG. 7: (Color online) The chemical potential that minimizes the binding energy for a prescribed charge in units 
of the fermion mass. 



TABLE I: The width, wh, of the minimizing Higgs profile at a prescribed charge, Q for different values of the 
Yukawa coupling constant /. 

estimate for the mass of the string is M ~ 5m? L. Taking m = SOOGeV and the length of the 
string to be the radius of the sun, L = Rq ^ 7 x lO^m we find M « 2.3 x lO^kg = W~^^Mq, i.e. 
only a very tiny fraction of the mass of the sun. On the microscopic scale, a string as short as 
the Compton wave-length of the heavy fermion would carry about 30 bound fermions and have an 
energy of slightly less than 9TeV. 



We have extended our previous spectral approach to find the leading quantum corrections to 
the energy of a cosmic string in a slightly simplified version of the electroweak theory. In the 
limit of many internal degrees of freedom, Nc — )• oo, these leading corrections come from fermions 
coupling to the string background. In this scenario Nq merely appears a combinatoric factor, which 
is justified by the asymptotic freedom of QCD. We have shown how to compute the distortion of the 
Dirac spectrum in the string background, and how to extract the full non-perturbative renormalized 
vacuum polarization using perturbative counterterms and conventional renormalization schemes. 
Substantial refinements of our previous techniques were necessary to make this calculation feasible, 
and we have presented a complete account including technical details in the appendices. Though 
we have focused on the computational method underlying previously published results, we have 
also discussed some novel results concerning the structure of the stable configuration. 

The basic idea of the quantum stabilization of cosmic strings is that the appearance of (near) zero 
modes in the distorted Dirac spectrum could help to produce negative contributions to the energy 
that overcomes the classical energy necessary to form the string. We have shown, however, that 
the contribution from the distortion of the remaining parts of the spectrum, i.e. the scattering 
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states, neutralizes the binding effect of the low-lying modes, resulting in a very small vacuum 
polarization energy. In particular, the vacuum is stable against spontaneous formation of weak 
strings for parameters that are physically sensible. The situation is more favorable for charged 
strings with explicitly occupied bound states, since such configuration need only be lighter than 
the same number of free fermions, to be stable on time scales over which we may neglect fermion 
number nonconservation. This approximation is valid if the inter-generation quark mixing is tiny 
(we have assumed it to be zero), and if the quark masses in the heavy fermion doublet are nearly 
degenerate (we have assumed exact degeneracy). 

For otherwise realistic parameters, we have shown that this binding mechanism sets in at sur- 
prisingly low fermion masses of around 300 GeV. This corresponds to values of the Yukawa coupling 
that are still small enough for our calculations based on the Standard Model to be reliable. A stable 
charged string can thus be formed when enough charge density of a heavy fermion doublet with 
about twice the top quark mass is available. 

If taken at face value, our findings suggest that a weakly coupled fourth generation of heavy 
quarks would make its footprint through the electroweak string phenomenology mentioned in the 
introduction — or conversely, that the non-observation of electroweak strings would put severe 
bounds on the masses of possible heavy quarks. However, such conclusions must be qualified by a 
number of simplifications that were necessary to make the calculation feasible. Most notably, the 
restriction to fermionic quantum fluctuations, although justified by the large A^'c-argument, leads 
to a quantum theory that is not asymptotically free, and in turn to the Landau pole problem at 
small string widths. We have presented a crude way of estimating this contribution in order to 
ensure that our findings are not affected by it. This treatment should obviously be improved by 
a full quantum calculation of the bosonic contribution to the vacuum polarization energy. Recent 
studies [231] in models related to ours indicate that the bosonic contribution can give interesting 
and non-trivial effects. We are currently investigating such an extension of our model. 

Other shortcomings of our model are the mass degeneracy of the fermion doublet, the lack of 
inter-generation quark couplings and, in particular, the decoupling of U{1) hypercharge. While the 
stable configuration that we have constructed can be embedded in a full SU{2) x U{\) model with 
multiple generations (since additional degrees of freedom would only serve to lower the energy 
in our variational approach), it is unclear if the new couplings provide new decay channels, in 
particular when bosonic fiuctuations are taken into account. We also plan to investigate such a 
scenario. 

The charges along the string may carry currents 39l44ll|. which in turn can have interesting 
consequences for baryogenesis and cosmology 24, |42]. This situation is similar to the Witten 
model [43|, |4J] and its generalizations, where the currents are induced by the coupling of extra 
scalar fields to the vortex. In this scenario, the Brownian network of vortices produced in an 
earlier (GUT-scale) phase transition contracts as the universe cools down. This process could 
eventually be stopped by the currents becoming superconducting, with the irregular vortex shapes 
being smoothed out by the surrounding thermal background to form circular rings. The final 
evolution stage would then be a universe filled with microscopic superconducting, charged vortex 
loops. Such a vorton [45|, |46|] universe has recently attracted much attention because it provides a 
viable candidate for dark matter with rather accurately computable properties that put stringent 
restrictions on cosmological models. It would be very interesting to study such a possibility in 
the electroweak standard model, with currents produced directly from fermions (as they are in 
our calculation), rather than from extra scalar fields. Although our present investigation does not 
directly address this question, it seems conceivable that a stable vorton could be created without 
requiring exceedingly large couplings or unrealistic masses. Combining this scenario to our picture 
could be another avenue for future research. 

Finally, it would of course also be interesting to study the Brownian network of strings as it is 
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TABLE II: Matrix elements of i<j ■ p. The 
table should be read as ( — h | icr - p | + + ) = 
1, for instance. The matrix is hermitian in 
combination with the derivative do- 
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TABLE V: Matrix elements of (cr • ^) (A) . 
The subscript denotes the argument of the 
trigonometric functions. 



formed in the phase transition if enough fermion charge is available. Due to their complexity, such 
configurations must presumably be studied in an effective (lattice) model. The necessary string 
interactions could potentially be addressed through further extensions of the spectral method. 
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Appendix A: Eigenvalue Problem 



To find the bound state spectrum, states with energy eigenvalues |e| < m, we first diagonalize the 
Hamiltonian matrix in the absence of a background potential, and then use these free eigenstates 
(with the proper boundary conditions built in) as a basis in which to compute the matrix elements 
of the background potential induced by the string. The diagonalization of this full Hamiltonian 
matrix in turn yields the fully interacting bound state spectrum. Note that in this procedure all 
states appear as "bound" states, since the volume of the coordinate space is finite. In general the 
energy eigenvalues of such "bound" states depend on the volume. However, the true bound states 
with |e| < m do not show finite size effects if the volume is chosen large enough because their 
wave-functions are located in a small sub- volume. 

The single particle Dirac Hamiltonian couples spin (5) and weak isospin (/) degrees of freedom. 
We can combine these degrees of freedom by introducing the grand spin states given in eq. ()12p . 
The matrix elements of the (two-component) operators entering the Hamiltonians in eqs. (jlOp and 
(|38p are listed in tables HIl to IVIII In all of these tables, we use the abbreviation Sx = sin(x) and 
Cx = cos(x), where the arguments of these trigonometric functions appear as subscripts. 
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TABLE VIL Matrix elements of Hp. The 
extra factor i leads to anti-hermitian matrix 
elements, compensating the same property of 
its spinor coefficient Pj^. 



Next, radial functions are introduced via the four-component spinors, cf. eq. ()15p . 



Je+n{kp)\l + +) 
KJe+n+i{kp)\i - +) 



{p\ei + 
{p\ei- 
{p\e£- 



Mkp)\e + - 

Kji+i{kp)\e - 
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Je+i{kp)\£-- 

K,Ji{kp)\e + - 
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U(p)i^ 




(h{p)\i 


- +) 


V53(P)|^ 


+ +) 
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WpW 


+ -) 



(Al) 



We note that k = sgn(e) y jip^ is weU defined for either sign of the energy eigenvalue since |e| > m. 

Using the dispersion relation for real momenta = k"^ + m^, we may also write k = = ^-jp- 
These two expressions are odd in k and thus suitable for analytic continuation k — )• it. The spinors 
to the left of the arrows in eq. (|Aip involve ordinary Bessel functions which solve the free Dirac 
equation. They will be used to construct the basis states for the Hamiltonian matrix. 

In the free case, the four spinors in eq. (|Aip each solve the Dirac equation individually, i.e. 
they do not couple. Once the background potential from the string is included, however, the radial 
functions fi and gi become distorted and mix under the dynamics. 

We now construct a discrete basis built from solutions of the free Dirac equation. To this 
end we must first impose the boundary condition that no flux runs from the center of the string 
through a circle at a large distance R. Since the flux is bilinear in the spinors with all products 
involving both an upper and a lower component, the no-flux boundary condition is equivalent to 
the requirement that either component vanishes. For the string winding n = 1, this amounts to 
the simple statement 



0. 



(A2) 



This conditions selects discrete momenta kf- in each angular momentum channel i, where r = 
1,2,... enumerates the momenta and thus the free basis states. We note in passing that a string 
winding n > 2 would require two separate sets of discrete momenta. 

The normalization of the spinors can be worked out using the Bessel function identity 

"1 



(A3) 



where A^*^^ are the roots of the Bessel function J^. Using furthermore the recursion relations for 
Bessel functions and their derivatives, we arrive at the following explicit expressions for the (free) 
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radial functions in equation (|T5|) . 



fl"{p) = N}"Je+i{K.p) 
f!i\p) = Nf\j,{krp) 



gPip) = NrJe+2ikrP) 

g!^\p) = NPje+,{krp) 

gt\p) = Nl'^J,+^{krp) 

g^P{p) = N^'\jeikrp) 



(r) 



ft\p) 



Nfh,+2{krP) 



(A4) 



where the superscripts on the momenta are omitted. The normalization factors are given explicitly 
by 



N 



(r) 



1 



1 



R\Ji+2{krR)\ 



e,. + m 
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1 sgn(er-) 
R\Je+2{krR)\ 



m 



(A5) 



To limit the number of basis states, we introduce a cutoff A and only include momenta with kr < A. 
This defines rmax, the maximal number of discrete momenta, which depends on R for fixed A. Due 
to the energy degeneracy this cutoff truncates the label r on the energy eigenvalues to run from 

1 . . . 2rmax) 



^2 



1 



max 



k?:_„ + 

' ' iimax 



(A6) 



+ l,...,2r„ 



Putting all the pieces together, we can now present the full Hamiltonian matrix i.e. the operator 
in eq. (llOp sandwiched between the spinors constructed above. The equations become simpler if 
we set 



C± = CAdG + Cgttg i ttr 



and 



S = SACtG ~ s^a^ 



(A7) 



with the right hand sides containing the elements of Vi in eq. (|17p . They are specified in terms 
of the string profile functions in eq. (|B4p of the following appendix. The interaction Hamiltonian 
matrix elements read 



(Ir 


Hint 


Is) = 


(2r 


Hint 


Is) = 


(3r 


Hint 


Is) = 


(4r 


Hint 


Is) = 


(Ir 


Hint 


2s) = 


(2r 


Hint 


2s) = 


(3r 


Hint 


2s) = 


(4r 


Hint 


2s) = 


(Ir 


Hint 


3s) = 


(2r 


Hint 


3s) = 


(3r 


Hint 


3s) = 



-at' 9? 



S(fi'^g['^+gVft'' 



-C+fi'-^g['^ 
S{ft'f['^+g^'g? 



ap fl 'g{ 



gt'f['' 



(r) [s] 
92 92 



+ S[i'\^2^+9Pf? 



ri Ar) A^) , n W {s) , ( Ar) (s) [r) As)\ 

C+fl /2 + C^gl 'g\ + ap l^f^ 'g\ - ffg V2 J 



S{ft'f?+9t'9'i' 



I Ar) (s) 

ap /2 53 



/ Ar) A^) 



(r) (s) 

93 9l 



27 



{3r\H,^,\As) = -C^ft\^^-C.gt^f^:'^ 

{Ar\H,^,\^s) = an (/f Vl^^ - gt^^^) + S [ft^gi'^ + ^f/i^^) . (A8) 

To keep the presentation simple we have omitted the radial integrals on the right hand sides, i.e. 
they are understood to be integrated with pdp (...)• total this defines Srmax x Sr^ax matrix 
elements of the interaction Hamiltonian. To populate the full Hamiltonian matrix, we set 

^ii2._,. = (2r|i/int|l5) (A9) 
for r, s = 1, . . . , 2rmax, and er+2qrn,^^ = Cr ioi q = 1,2,3. This yields the 8rmax x Sr^ax matrix 

Hr,s = €r6rs + H^^) r,S = l,..., Sr^ax , (AlO) 

which is diagonalized numerically by means of a Jacobi routine. 

Once the radius R and the momentum cutoff A are large enough the true bound state spectrum 
should become stable against further increase of these parameters. Typical values are A ~ 8m 
and R = 75/m, so that our free basis comprises about 400 energy eigenvalues, each with fourfold 
degeneracy and the Hamiltonian matrix has 1600 x 1600 entries for the lowest angular momentum. 
For wider string profiles, bound states occur in higher and higher angular momentum channels. For 
instance, bound states appear for up to ^ = 5 when wh ~ 6/m, while narrow widths wh < 1/m 
only induce bound states in the channel £ = —n = —1, i.e. the effective S-wave channel. 

We have verified the gauge independence of the bound state energies by checking that they 
are insensitive to variations in the shape of the gauge transformation profile ^(p). Also the zero 
mode in the £ = —n = — 1 channel is observed for = '^/2 regardless of the values of the width 
parameters. 



Appendix B: Scattering problem 

In this appendix we describe the scattering solutions to the Dirac equation (fT7|) . To this end we 
write the Dirac Hamiltonian, eq. (jlOp in terms of 4 x 4 matrices and derive the differential equation 
for the Jost function. 

1. Differential equation for Jost function 

The derivative operators as well as the angular barriers are contained in the diagonal 4x4 
matrices 

Du = dpi-- diag( -{i + n + l),-{i + l),i + n,i) 
P 

Dd = dpl-^dia.g{i + n, i, -{i + n + 1), -(^ + 1)) (Bl) 

where C = diag(— 1, —1, 1, 1). For the interactions, we must compute the matrix elements of the 
various pieces in eq. (fTOj) within the grand spin basis eq. p^ . The explicit expressions for the 
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emerging radial functions are listed in eqs. (jlSp and (jl6p . cf. eq. ()Aip . After some lengthy algebra, 
the interaction matrices in the system (I17p can written in terms of simpler sub-matrices, 



(B2) 



P \-P G. 



where the 2x2 submatrices are 



_ / sinA cos A \ / — sin^ cosM _^ / — 1 \ . . 

^ ~ ycosA -sin Ay ^^^y cose sin^y 1^1 y ' ^ ' 

The coefficients an, Op, qg, and ar are radial functions determined by the background profiles 
fci fn and the gauge function ^, 

1 dap) 



ar{p) 



2 dp 



71 

o^cip) = Trfcip) sinA(p) 
2p 

Tl 

a«(p) = ^ (/g(p) - 1) sine(p) 
anip) = m {fnip) cos A(p) - l) 

ap{p) = mfnip) smA{p) . (B4) 

Note that the new gauge function also enters via A(p) = — (,{p)- We now write the Dirac 
equation for the matrix fields defined in eq. (I25p as 



dpT = [Mff + Od]-J' + J'-M^J'} + [Mfg + kG]-g-Zd (B5) 

dpG = [Mgg + Ou] -g + g- A^g + [Mgf - kG] ■ F ■ , (B6) 
where the 4x4 matrices without an overline are purely kinematic, 

\He+n+i{kp) ' Hi>+i{kp) ' Hi+n{kp) ' / 

Zd = diag f ^i+n+i{kp) H(,+i{kp) ^ 
V Hi+n{kp) ' -ff^M ' He+n+i{kp) ' Hi+i{kp) J 

Ou = -diag (-(^ + n + !),-(£ + !),£ + n,£) 

Od = -diag (£ + n,£,-(£ + n + !),-(£ + 1)) 
P 

C = diag(-l, -1,1,1). (B7) 
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The matrices multiplying T and Q from the right are also independent of the background potential, 



M 



(r) 
// 



-kC -Zd-Od and = fcC • Z„ - 0„ . 



(B8) 



Genuine interactions from the string background are solely contained in the overlined matrices in 
eq. ()B6p . Using the same 2x2 matrix notation as above, we have explicitly 



Mgg = CVud 



G+ P 
-P -G. 
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ff 



-CV„ 



du 



M 



gf 



„ ^ y 11.11. — „ 
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-nCVdd = K 



-G_ P \ 
-P G+) 

-H G_ 



(B9) 



-G+ H 



The solutions to the differential equations (|B6P subject to the boundary conditions T ^ 1 and 
^ — >■ ][ at /3 — 7- oo define the scattering solution, eq. (j27p . from which we extract the scattering 
matrix as described in eq. ([29]) . 



2. Born series 



To set up the Born series defined in eq. (|30j) . we simply expand the system of differential 
equations from the last section in powers of the background potential, which only enters the 
overlined matrices. At first order, we obtain 



dpT^^^ = Od ■ -F(i) + ^(1) • M^pj + kC ■ -Zd + Mff + Mfg-Zd (BIO) 

dpG^^^ = On ■ G^^^ + g^'^ ■ M^pg - kC ■ J-(i) • Z^ + Mgg + Mgf ■ Z^ . (Bll) 

The matrices M^^^ do not contain the interactions and are thus of order zero. In the same way we 
obtain the second order equations, 

= On ■ Q^"-^ + a^') • Mp;} - kC ■ . Zu + Mgg ■ Q^^^ + Mgf ■ -F(i) • Zu . (B13) 
With these Jost-like matrices, the Born series for the 5-matrix is 5 = 1 + S^^^ + S^'^^ + ... with 

5« = lim{7^-i.[j-«*-J-«].K} 

5(2) = hm {n~' ■ [j-(i) • (j-(i) - J-(i)*) + - ^(2)] . n:} (B14) 

and similarly for Gi with Ti^ — s- Tid- The Born expanded eigenphase shifts are now simply given, 
to first and second order, by 



.(1) 



1 



tr 



Im [SI 



and 



?(2) 



■tr 



(B15) 



The third and fourth order pieces will be treated as part of the fake boson formalism discussed 
below in section fC 2[ 
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3. Analytic continuation 



We describe the continuation to imaginary momenta for the case where e = + ni^] the 
second Riemann sheet (e = —s/k"^ + m?) works analogously. The analytic continuation concerns 
the Hankel functions, which turn into modified Bessel functions, and — t- Yd, with 

V ( Kg+njtp) Kjjtp) Kt,+n+i{tp) Ke+i{tp) \ _ _ .-i ^ 

Yu - diag , , KeUtp) ' Keitp) J " ^"""^ ' ^""'^^ 

Furthermore, the kinematic coefficient turns into a pure phase 



771 + i\/t'^ — m? 

t 

The system of differential equations for imaginary momentum then becomes 



(B17) 



dpF = [Mff + Od]-F + F-Mff + [Mfg-tC]-g-Yd (B18) 

dpQ = \Mgg + Ou]-Q+Q- M^jJ + [Mg f + tC] ■ F ■ Y^ . (B19) 

with the boundary conditions that F and Q both approach unity at /) — )■ oo. For simplicity, we 
have omitted the momentum arguments in the radial wave-functions J- and G and also used the 
same symbol as in the case of real momenta, eqs. (|B6p . The coefficient matrices in the differential 
equations are slightly modified: 



M^J'J = -tC-Yu-Ou Mfj =tC-Y, 





(B20) 



while Mgg and M/f are the same as on the real axis. 

Unlike the Schrodinger problem, the differential equations in the present case do not become 
real on the imaginary axis. Rather, charge conjugation e — t- — e induces complex conjugation. It 
is therefore not surprising that the naive extrapolation limp_;,o det(J^) does not give a real result. 
Instead, we find numerically that F = G* with the imaginary part being independent of angular 
momentum for a given value of t. The origin for this imaginary part lies in the subtle definition of 
the Jost function via the Wronskian between the Jost solution, i.e. F or Q, and the regular solu- 
tion, which satisfies momentum independent boundary conditions at the origin. The momentum 
independence of these boundary conditions ensures that the regular solution is an analytic function 
of complex momentum. Analyticity of the Jost solution, on the other hand, is guaranteed by the 
non-singular behavior of the interaction potentials, which is in turn a consequence of the bound- 
ary conditions on the profile function £,{p). At the origin, the Higgs field differs from its vacuum 
expectation value (it actually vanishes), which modifies the relative weight of the upper and lower 
Dirac components. More precisely, the non-diagonal elements of the matrices in eq. ()B20p vanish 
at the origin and the eight differential equations decouple with respect to the spin and weak isospin 
index on the radial functions in eq. (I16p . For real momenta k, a typical solution in the vicinity of 



p = then looks like [22] 




( ^/E + mcAfH{0)Ji{qp) \ 



(B21) 



- mcAfH{0) Ji+i{qp) 
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with q = \JE'^ — (mc/^fHiO))'^ and similar dependencies for the other six radial functions. The 
square-root coefficients cause the proper definition of the logarithmic Jost function, to be 

exp [iy{t)] = ( , \ lim det(^) = ( ^ ^ , lim det(g) (B22) 

^^^'^ \T - imcAfHiO) J P~^o ^ ' \T + imcAfHiO) J P^o ^"'^ ^ ^ 

with T = \/t'^ — m?. The power of two occurs because we compute the determinant of a 4 x 4 
matrix. Notice that this redefinition not only cancels the imaginary parts, but also modifies the 
real part. Furthermore it avoids the logarithmic singularity in In [limp^o det(J^)] otherwise observed 
numerically at t ~ m. Since fu is part of the interaction, the correction prefactor in eq. ()B22p also 
contributes to the Born series. To make this explicit, we write 

\ ( T — im \ T — im \ fB23) 



ianiO) 1 / «_f/(0) ^ ^ 



im 2 \T — im 



+ . . . , (B24) 



and subsequently set 0^(0) = —m. The Born expansion of the remaining determinant in eq. ()B22p 
is constructed as for real momenta by iterating the differential equation (IB19P in the interaction 
Mi. 

Numerically, we integrate the differential equations (|B6p . ()B19p . their Born expansions and the 
fake boson analog^ 



Jpe{t, p) = 2tL,{tp)dpV,{t, p) - ujit, p) + V{p) with L,{z) = " ^ 



from some large radius /5max 4po to pmin ~ with the boundary condition T{praax,f^) = 
and identify limp^Qj^(p,k) = T{p^\^,k). Alternatively, this identification can also be obtained 
from the derivative of the wave-function. Furthermore, a differential equation is formulated for 
lndetJ^(/3, A;)* to avoid 2tt ambiguities in the computation of the phase shift, (5^(A;), c/. 
eq. (IMD. 

The computations for real momenta have been performed mainly for use in the consistency tests 



on the unitarity of the scattering matrix and the spectral sum rules SlJ. There is one more merit of 
considering real momenta: Channels that include Hankel functions with zero index [I = —2, —1, 0) 
are particularly cumbersome because regular and irregular solutions are difficult to separate in 
such cases, because they go as a constant and ln(/3), respectively, at pmin ^1, c/. eq. ([29]) . As a 
consequence, pmin niust be taken tiny in the problematic channels to obtain the correct scattering 
matrix in eq. (j29p . On the real axis, the result can be checked against extracting the 5-matrix 
from the derivative of the scattering wave-function because Yq{p) ~ diverges like a power. 

For calculations on the imaginary axis, we assume p^mn ^ 10 and successively carry out an 
extrapolation 

_ai ^ 02 

In(pmm) ln^(p n 

for the Jost function in these channels. We test the final result, i.e. vq, for stability against 
further changes of pram and also check the condition Im(fo) = 0. In the non-problematic channels 
I ^ {—2, —1, 0}, it is sufficient to set pmin ~ 10^^^ in order to represent the origin. 

We have also successfully tested our numerical results of the scattering data against the reflection 
symmetry £ — t- — + 2n). 



I'iPmin) = J^o + 7777 \ + T~27" T ■ ■ ■ ' {B2Q) 



The boundary condition is v tit, oo) = dpUt{t,oo) = 0. The second order contribution required in eq. 
obtained from the expansion Vi = + ^^f^' + ■ • ■, where the superscript labels the order in V{p). 



32 



Appendix C: Feynman diagrams 

In this appendix we describe the details of the computation of the Feynman diagrams. We start 
from the series in equation (j4ip . This computation involves three parts: 

1. The contribution linear in Hj. This term vanishes identically with the no-tadpole condition. 

2. The piece quadratic in Hj. This term is quadratically divergent at high momenta and must 
be carefully regularized to handle the leading and subleading logarithmic divergences. 

3. The contribution from terms cubic and quartic in Hj. They are only logarithmicly divergent 
which makes the separation of the finite parts simpler. Since the corresponding Feynman 
diagrams are complicated to evaluate we employ the fake boson methods to compute this 
part of the vacuum polarization energy. 

We first consider the MS scheme in which only the bare divergences proportional to 




are subtracted, and then determine the finite counterterm coefficients suitable to implement the 
on-shell renormalization scheme. 



1. Second order contribution 



After imposing the no-tadpole condition^ we find the contribution to the action functional up 
to second order in Hj within the MS scheme as 



1 

"8^ 



d'^k 



dx In 



(2vr)4 

xtr/|[m^ - x{l - x)k'^^ 



1 — x{l — x) 
-L{k) ■ L{-k) - 3 {h{k)h{-k) +p{k)p{-k)) 



+x(l — x) 



k ■ L{k) k ■ L{-k) - -k^L{k) ■ L{-k) 



+- 



+2rn?p{k)p{—k) + imk ■ L{k)p{—k)^ 
"^^^ in [h{k)h{-k) + p{k)p{-k)\ . 



87r2 J (27r)4 6 



(CI) 



Here the fields with momentum arguments are the Fourier transforms of the corresponding spatial 
fields in eq. (pQl) . Specifically, we introduce the notation k := k^ = (A:o, A;_Lk^ + k^z)^ and L := 
L'^ = (Lq, L_l + Lsz)^ with Lq = = 0. As a result, we have 



h{k) 
p{k) 



h(-k) = -i27Tf5{ko)6{k;)ho{k^) 

p\-k) = -{27rf6{ko)6{k3)i-irpn{k^)Ip{'Pk) 



(C2) 



^ The C3 type counterterm in eq. (|43p contains a term quadratic in the fluctuations about the Higgs vev. Its finite 
contribution is essential to keep the pseudo-scalar part of the Higgs field massless, i.e. the expansion of the 
coefficient of p{k)p{ — k) starts at 0{k^). 
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i=l 

where (pk is the azimuthal angle in momentum space. The coefficients are isospin matrices, 



(C3) 



l^^ = 



if = {-ir-^a\^'{k^)Ip{^k) 



(2) 



iasikj 



-1 0' 
1, 



if = i-ira'-'{k^)Ip{ipk) 



-1 0' 
1, 



(C4) 



The matrix Ip is defined in eq. (I40p . here to be taken as a function of the azimuthal angle in 
momentum space. The functions Hq, p„, a 



(-) 



are the Fourier transforms 



hoik) 

Pn{k) 



pdp aH{p)Jo{kp) 



pdpap{p)Jn{kp) 



a^^\k) = / pdpar{p)[Jn+i{kp) ± Jn-i{kp)\ 





as{k) 
af\k) 



pdp [aG{p)s/^{p) - ag(p)se(p)] Jii^p) 
pdp [aG(p)cA(p) + a(^{p)ci:{p)] [Jn+i{kp) ± J„_i(/cp)] 



Some of these terms can be conveniently combined, 

3 

Ip{^k) 



i=l 
3 



n-1 



E4 



-1 



-1 0' 
1, 



Then we find the second order contribution to the energy 



(2) 
FD 



kdk k ^ ^ 4^2 j^p2 ^ 2mA;/i (q(+) + a^'A Pn 



Att 3 



+k^h 



(m^/i + k^I^) 



(as)' 

2 



with the Feynman-parameter integrals (r/ = /c/m) 

/■I 2 

h = dxlnll + x{l - x)r]'^] = -y% + rf arsmhCnh) - 2 , 
Jo V 

/I ^ /4 _|_ jn2 2 

- x) In [l + x{l - x)r]'^] = ^ ^ 3 [r]"^ - 2] arsinh {v/2) + 



(C5) 

(C6) 
(C7) 



(C8) 



3r/2 18 



• (C9) 



34 



2. Fake boson method 

We have already discussed the spectral part of the fake boson approach in eq. (IB25P . Here we 
focus on the Feynman diagram part. First we need to determine the logarithmicly divergent con- 
tribution to action from the third and fourth order Feynman diagrams. They can be parameterized 
by a radial integral, 

cp = J pdp ^[ajj + a^,) [ajj + aj:) + Aman) + 'ictr {ctHOi'p — oipa'fj) 

+4 (a^ + Oq + a| + la^acc^j^) {ajj + Op + Iman) 
64 2 / 2 2 \ 8 „, 

-4-ap [an (ogca + a^c^) + ap (agSA - o:^s^)] | , (CIO) 

where primes denote derivatives with respect to the radial coordinate. With this radial integral 
the divergence reads, in dimensional regularization. 



m 



(Cll) 



Here T and L are the (infinite) lengths of the time and z-axis intervals, respectively. 

A boson field that fluctuates about a background potential V{p) = nn?-^^e~'^P/ P'-^ causes a 
similar logarithmic divergence for its vacuum polarization energy at quadratic order. In fact, the 
only replacement in eq. (IC11|) is cf —)■ c_b with 

cs = \ r pdpV\p) = ^-^. (C12) 



4 7o 512 

As for the spectral part, eq. (j48p . we rescale the fake boson potential with the strength of the 
fermionic divergence cp so that we are only left with the finite part of the second order (boson) 
Feynman diagram, 

cp /■°° kdk ^ ^^2 



AEb = -— tttW- (C13) 
CB Jo 167r 

In this equation, the Fourier transform of the fake boson background is 

Voik) = pdpV{p)Jo{kp) = m'pl (C14) 

■^0 [4 + k^pl] 2 

For the numerical test mentioned after eq. (j50|) we vary po and verify that the vacuum polarization 
energy does not change. 

3. On shell renormaUzation 

We parameterize the counterterm coefficients in dimensional regularization 

— I / T^--^(/2-l + k)"^ + c,, (C15) 
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for s = 1, . . . ,4. In the MS scheme, ci = C2 = C4 = 0, while the on-sheh conditions discussed in 
the main text yield 



ci 



C2 



C4 



12 (47r)2 
2/2 



21n [1 — x(l — 



2(47r 



1 + 6 / dxx (1 — x) 
^0 

|- + 6y dxx (1 — x) In [1 — — x)^^] I 
|;uf^+6 y dx In [1 — x(l — 



w\ 



1 — x{l — x)iJ? 



w 



(C16) 



We recall that the no-tadpole condition implies C3 = / {A-k"^ v^) = f'^m? / {A.-k'^). In this scheme, 
the pole position of the gauge boson field is not prescribed but rather becomes a prediction. We 
find an implicit equation for the gauge boson mass = my/ /m: 



2/2 ^ 16^2 I 3 



x2(l -x)2 

1 - x(l - x)//^ 



-6^ dxx{l - x)ln [1 - x{l - x)nly] - dxln [l - x{l - x)^^] | . (C17) 



Appendix D: Landau ghost estimate 



In the present treatment (without gauge boson loops) our model is not asymptotically free. This 
results in unphysical poles of the renormalized propagators at large space-like momenta. These 
so-called Landau poles are not real singularities but rather indicate the breakdown of our treatment 
in certain momentum or parameter regimes. In the present model the problem has a notable effect 
only for narrow background profiles and/or large coupling constants. We have implemented a 
procedure similar to that of ref. [47| to verify a posteriori that the interesting configurations do 
not suffer from this unphysical effect. 

Specifically, we write the renormalized quadratic contribution to the energy per unit length 
coming from the pseudoscalar component of the Higgs as 



_fq_ 

(2vr)2 



tr 



p{q)p{-q) G~\q^) 



which involves the corresponding (inverse) propagator for space-like momenta. 



2 ^ f'Nc 



dx x(l — x) In 



m2 + x(l — x)(/2 
m? — x(l — x)m'\j 



-2m2 / dxx{l—x)\u. 
Jo 



1 + x(l - x)- 



(Dl) 



In the vicinity of the Landau pole (g2 ~ ttiq) this propagator has the expansion 



(D2) 
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where 



dip- 



(D3) 



is the residue of the pole. This allows us to remove the Landau pole explicitly by introducing 



Zg 



m 



G 



(D4) 



We eliminate the artificial ghost contribution associated with the Higgs field from the energy in 
chirally symmetric way 



E. 



G 



(fq 



1 



(D5) 



To study the effect of the Landau ghost, this quantity should be compared to the same contribution 
without the Landau ghost removal in eq. ()D4p . which we call E^^\ 

In the same way, we can treat the gauge boson contribution to the renormalized energy per unit 
length, 



_fq_ 

(2vr)2 



tr 



W^Aq)W 



(D6) 



where W^i,{q) denotes the Fourier transform of the field strength tensor for the static background, 
eqs. ([5]) and dH]), while 



G^{q^) = l + 



Nc9^ 
167r2 







^ +x{l- x)q'^ 

ax x[ \ — X) In — - 



x(l — x)m 



2 

w 



(D7) 



describes the inverse gauge field propagator for space-like momenta. Again, this propagator has a 
pole at q^ = tuq with residue Zq which we remove by defining the subtracted inverse propagator 



1 



G^'{q' 



Zg 

- 2 - 2 



(D8) 



The Landau ghost eliminated gauge field energy then becomes 



E. 



{W) 
G 



(2vr)- 



■tr 



W^.MW^'i-q) AT^Hq 



(D9) 



Asymptotic freedom implies that the Landau poles at large spacelike momentum in the various 
propagators should disappear at any order in perturbation theory, and also in the full theory. We 
therefore expect that the difference between E^^^ + and E^^ + Eq^^ is small whenever the 
effect of the unphysical Landau ghost in our model can be safely ignored. For the model parameters 
that we found interesting, g = 0.72 and / ~ 2, this condition is indeed satisfied since the relative 
difference between E^^'^ + and E^^^ + E^^^ is only a fraction of a percent even for narrow 

configurations with wq = wh = 0.1. 

To see that there are indeed background potentials where the Landau ghost contribution is 
sizeable, we present the same comparison between E'^^^ + s(^) and E^^"^ + E^^^ for g = f = IQ 
in table IVIIII We observe that the Landau ghost causes the well-known instability for narrow 
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WG WH 




+ 4^) 


0.1 


-15.597 


3.220 


0.5 


-0.168 


0.209 


2.0 


0.041 


0.082 


4.0 


0.061 


0.077 


6.0 


0.068 


0.077 


8.0 


0.070 


0.077 



TABLE VIII: Landau ghost removal for g = f = 10. As an example we have chosen = O.Stt. 



configurations and large couplings |48l ]. However, for wider configurations its effect is moderate 
even when the coupling is large. It should be emphasized that the present approach to the Landau 
ghost problem is only qualitative since the energy expressions ()D2p and ()D6P are not rigorous. 
However, the present method convinces us that the configurations discussed in the main body of 
this article do not suffer from this problem. 
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